Phys. Rev. A 82, 053606 (2010) 



Three-dimensional gap solitons in Bose-Einstein condensates supported by one-dimensional 

optical lattices 

A. Munoz Mateo 1 , V. Delgado 1 , and Boris A. Malomed 2 

1 Departamento de Fisica Fundamental II, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain and 
2 Department of Physical Electronics, School of Electrical Engineering, 
Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel 
(Dated: 11 June 2010) 

We study fundamental and compound gap solitons (GSs) of matter waves in one-dimensional (ID) 
optical lattices (OLs) in a three-dimensional (3D) weak-radial-confinement regime, which corresponds 
to realistic experimental conditions in Bose-Einstein condensates (BECs). In this regime GSs exhibit 
nontrivial radial structures. Associated with each 3D linear spectral band exists a family of fundamental 
gap solitons that share a similar transverse structure with the Bloch waves of the corresponding linear band. 
GSs with embedded vorticity m may exist inside bands corresponding to other values of m. Stable GSs, 
both fundamental and compound ones (including vortex solitons), are those which originate from the bands 
with lowest axial and radial quantum numbers. These findings suggest a scenario for the experimental 
generation of robust GSs in 3D settings. 
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I. INTRODUCTION 

A ubiquitous tool for the control of collective excitations 
in Bose-Einstein condensates (BECs) is provided by opti- 
cal lattices (OLs), which are induced by the interference of 
laser beams illuminating the condensate [1]. OLs are espe- 
cially efficient in supporting matter- wave solitons. It has 
been predicted that two- and three-dimensional (2D and 
3D) OLs can stabilize solitons against collapse [2]. OLs 
acting in the combination with repulsive interactions can 
give rise to diverse species of gap solitons (GSs), in one di- 
mensional (ID) y,|4[] and multidimensional [5, 6] geome- 
tries. Quasistable GSs were created in the 8T Rb condensate 
loaded into a cigar-shaped trap incorporating an axial OL 
[7]. Extended states, built as segments of nonlinear Bloch 
waves trapped in the OL, were also reported [8]. 

Most theoretical studies of GSs have been carried out 
in the quasi- ID regime, assuming that the transverse con- 
finement is tight enough to reduce the wave function in 
the transverse plane to the ground state of the correspond- 
ing 2D harmonic oscillator (HO) fid. 

In this case, 

the description of the relevant dynamics amounts to the 
ID Gross-Pitaevskii equation (GPE) in the axial direction, 
which follows from the underlying 3D equation after aver- 
aging out the radial (transverse) degree of freedom. Since 
in the quasi- ID regime the radial quantum huu± is the 
largest energy scale of the system, GSs cannot decay by 
exciting higher-order radial modes, being therefore partic- 
ularly stable. However, the creation of genuine quasi- ID 
settings, strongly squeezed in the radial direction — for in- 
stance, those realizing the Tonks-Girardeau gas [11] — is 
a challenging experimental problem. The settings used in 
experiments with GSs 01, 0, S] actually correspond to the 
weak transverse confinement (see details in Sec. IV). 

In this work we aim to predict matter- wave GSs in the 
regime of the weak radial confinement, characterized by 



a recoil energy of the axial OL, E R , comparable to hoj±. 
While this regime is most relevant to the experiment, the 
formation of GSs under these conditions was not yet stud- 
ied theoretically. For instance, E R /huu± = 1 corresponds 
to the 87 Rb condensate, with s-wave scattering length a s — 
5.29 nm, confined by the combination of the transverse 
trapping frequency (jJ±/2tv = 240 Hz and axial OL of pe- 
riod d = 1.55 /im (physical results given in this article 
correspond to this typical setting). In this regime, the axial 
GS structure may readily excite higher modes of the radial 
confinement; hence the 3D character of the dynamics is es- 
sential and the ID reduction cannot be used. The situation 
is somewhat similar to that for quasi- ID GSs, which were 
predicted, in the framework of the density-functional de- 
scription, in fermionic superfluids [12]. In that case, the 
underlying Fermi distribution implies the filling of many 
transverse energy levels. 



While there are models that generalize the ID GPE 
by taking into account small deviations from the one- 
dimensionality Jsl, [Toll , we consider the setting in which 
the axial and radial directions are equally important and 
inseparable; hence the use of the full 3D equation is nec- 
essary. We demonstrate that stable solitons, which are true 
gap modes in terms of the underlying 3D band-gap struc- 
ture, exist in this regime, suggesting possibilities for the 
creation of robust 3D solitons. This objective is of princi- 
pal significance because, thus far, no truly 2D or 3D matter- 
wave solitons, nor their counterparts in optical media with 
the Kerr nonlinearity, have been created, in spite of many 
theoretical predictions [13]. We also find solitons inside 
the bands, which may exist due to the difference in the az- 
imuthal index between the soliton and the band. 
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FIG. 2: (Color online) Isosurfaces of the atomic density, taken 
at 10% of the maximum value, for GSs corresponding to points 
F1-F8 in Fig. \Qb)- Panel (i) depicts a GS in a weaker lattice, 
with s = 5. The respective quantum-number sets, (n, m, n r ), are 
indicated in each panel. The length scale of the left part in each 
panel (reduced to 45% in comparison with the image on the right) 
is gauged by the OL period. 



FIG. 1: (Color online) The band-gap structure, produced by 
the linearized GPE, with equal axial and transverse energies, 
E r /Hujj_ — 1 (white areas are gaps). The normalized chemi- 
cal potential is shown vs the scaled lattice depth, s (a), and as 
a function of quasimomentum q for s = 15 (b). The right-hand 
panels in (b) indicate the location of the GSs displayed in this pa- 
per. Sets (n, m, n r ) represent the quantum numbers, and the bold 
red curves depict the band-gap structure in the ID model. 



II. THREE-DIMENSIONAL GAP SOLITONS 

The 3D GPE, with the radial-HO and axial potentials, 

Vj_(i*_l) = (M/2)cj 2 r 2 and KO) = sE R sin 2 (irz/d), 
is 

ihi/> t = [-(fi 2 /2M)V 2 + V ± (t ± ) + V z (z) + gN\t/j\ 2 ^, 

(1) 

where N and M are the atomic number and mass, g = 
47Th 2 a s /M, d and s are the period and depth of the OL, 
the recoil energy is E R — {nKf / (2Md 2 ), and the norm 
of the wave function is 1 01411 . 

Figure Q2a) displays the scaled shifted chemical poten- 
tial, Jjl = {n — hou±)/E R , of the noninteracting 3D con- 
densate (g = 0) with E R /hoj± = 1, as a function of 
the OL strength s. In the plot, bandgaps separate shaded 
Bloch bands. For the sake of comparison, the band-gap 
diagram obtained from the corresponding ID equation is 
also shown (bold red lines). Since the energy levels of 



the radial HO are E = {2n r + \m\ + l)hu±, where 
n r = 0, 1, 2, . . . and m = 0, ±1, ±2, ... are the radial 
and azimuthal quantum numbers, the excitation of trans- 
verse modes gives rise to a series of replicas of ID Bloch 
bands, shifted up in the energy by multiples of hujj_/E R . 
In what follows, we present most of our results for a strong 
OL with s = 15, which corresponds to the dashed ver- 
tical line in Fig. [Ha). The respective spectrum is shown 
in Fig. QIb), which displays fl as a function of quasimo- 
mentum q in the first Brillouin zone. The 3D bands are 
characterized by quantum-number sets (n, m, n r ), where 
n = 1, 2, 3, ... is the band index of the corresponding 
ID axial problem. As in Fig. [Ha), the superimposed bold 
red lines represent the results generated by the correspond- 
ing ID equation. As indicated in Fig. Qlb), the lowest 
band corresponds to (n, m, n r ) = (1, 0, 0). The next two 
bands, with numbers (1, ±1, 0), have equal chemical po- 
tentials, being replicas of the lowest band shifted upward 
by hu_\_/E R . Likewise, the next three bands, with num- 
bers (1, ±2, 0) and (1, 0, 1), are shifted by 2Hcj ± /E r and 
so on. Dots in the right-hand panels in Fig. [Tfb) indicate 
the location of the GSs that we consider in this paper, with 
the horizontal axes indicating the number of 87 Rb atoms in 
each nonlinear state. 

GS solutions have been obtained as numerical solutions 
of the stationary version of Eq. ([T]), using the Newton con- 
tinuation method with the Laguerre-Fourier functional ba- 
sis. In Fig. [2 we display a set of solutions for GSs which 
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FIG. 3: (Color online) The Jl(N) curves representing the four 
lowest-lying GS families. The inserted panels display the soli- 
tons at marked points. The field of view in each panel is 
4.6 fim x4.6 /xm, in terms of the 87 Rb condensate. 

correspond to points F1-F8 in Fig. QIb). Stable GSs corre- 
spond to n r — [panels (a)-(c) and (h)]. In addition, Fig. 
Hi) displays an example of a loosely bound but also stable 
GS supported by a weaker OL, with «s = 5. The param- 
eters are adjusted to N — 450 atoms [N = 350 in Fig. 
Eti)]; in Ref. 10], the GS had ~ 250 atoms. The right-hand 
image in each panel is a perspective image, while the top 
and bottom left plots are axial and lateral views. 

As shown in Fig. [3l GS families are represented by 
Jl(N) curves, each approaching a certain band at N — » 0, 
whose set (n, m, n r ) is used to label the families. As N 
increases, the size of the GSs increases too, keeping a char- 
acteristic soliton structure. The GS with the lowest chemi- 
cal potential, of type (1, 0, 0), which corresponds to point 
Fl in Fig. [Hb), is shown in Fig. Eta). GSs of this type 
are disk-shaped objects localized within a single cell of 
the axial OL. While the radial shape of this GS approaches 
the ground state of the corresponding HO at N — » 0, the 
contribution of radially excited states becomes more pro- 
nounced with the increase of N, making the GS a fully 3D 
mode. In particular, point Fl in Fig. \Hb), with Jl = 6.37, 
corresponds to fi = 7. 37 'huj± in physical units, which 
is much larger than the radial-excitation quantum, Huj±. 
The decomposition of this GS over the basis of radial HO 
modes yields a mixture of states n r = 0,1, and 2, with 
respective weights 76%, 22%, and 2% 

As per azimuthal index m, the GSs in the (1,1,0) and 
(1, 2, 0) families carry vorticities 1 and 2; see Figs. [2]and 
O (The vorticity may also be embedded into solitons in the 
case of the tight transverse confinement, with hou± ^> E R 
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FIG. 4: (Color online) Gap-soliton complexes corresponding to 
points C1-C4 in Fig. DIb). The underlying ID lattice potential is 
also shown. 

fToh . as well as into GSs supported by the OL in the 2D 
geometry [6].) The GSs corresponding to n r > 1 feature 
a complex radial structure, composed of many HO modes. 
Nevertheless, these solitons exhibit a set of zero-density 
rings reminiscent of the HO wave functions with the re- 
spective values of n r \ see Figs. Etd), Of) and[2tg) for 
n r = 1, 2, and 3. 

The higher-order (n = 2) band of the axial potential 
gives rise to GS families of types (2, ±1, 0). As seen in 
Fig. (2h), they exhibit two major peaks in the axial direc- 
tion with a node between them, all squeezed into a single 
lattice cell (in the ID GPE, solitons of this type are known 
as subfundamental GSs [4]). On the other hand, because 
band gaps emerging at large values of the chemical poten- 
tial are very narrow, no GS of type (2, 0, 0) exists here, as it 
is not able to place itself within the corresponding narrow 
gap. It is worth noting, too, that GSs of the (1, 0, 3) and 
(1, 1, 2) types are found inside the (2, 1, 0) and (2, 0, 0) 
Bloch bands, respectively [the former case corresponds to 
point F7 in Figs. [Tfb) and[2tg)], which is possible be- 
cause modes with different azimuthal numbers m do not 
mix. In fact, these may be understood as examples of "em- 
bedded solitons", whose chemical potential falls within lin- 
ear bands. While they were studied in detail in ID models 
fl5ll . no example of multidimensional embedded solitons 
has been reported previously. 

Fundamental gap solitons play an important role as ele- 
mentary building blocks of higher order localized nonlin- 
ear structures [16]. Some examples of the latter, gener- 
ated from the symmetric linear combination of three fun- 
damental solitons, are shown in Fig. |4| The gap solitons 
displayed in Figs. @Ja) and@Jb), which are composed of 
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FIG. 5: (Color online) (a) The stable evolution of a GS of the 
(1,0,0) type after the application of the perturbation, (b) The 
instability of the soliton of the (1, 0, 1) type, from Fig. Od). 

three identical (1, 0, 0) and (1, 1, 0) solitons, respectively, 
and correspond to points CI and C2 in Fig. Qlb), have ex- 
actly the same chemical potentials as their respective el- 
ementary constituents (jl — 6.37 and 6.79) but contain 
approximately three times as many particles (N — 1372). 
The compound gap solitons in Figs. Hfc) andH^d), which 
have the same particle content as the previous ones but are 
higher in the spectrum (points C3 and C4), are bound states 
of fundamental constituents of different types. The for- 
mer is composed of one (1, 0, 0) and two (1, 0, 1) solitons, 
while the latter contains one (1, 0, 1) and two (1, 0, 2) fun- 
damental solitons. 



III. STABILITY 

The stability of the GSs was tested by simulating their 
perturbed evolution. To this end, Eq. (Q} was solved by the 
Laguerre-Fourier pseudospectral method, using the third- 
order Adams-Bashforth time-marching scheme. To include 
all potentially dangerous disturbances, we perturbed the 
GSs by simultaneously increasing the OL depth and de- 
creasing its period, both by 2%, translating the OL in the 
axial direction by 2% of the lattice period, and applying a 
2% quadrupole deformation to the transverse trapping po- 
tential. After waiting for t — 1 ms, the combined pertur- 
bation was removed, allowing the system to evolve for 2 
s. The simulations demonstrate that the fundamental GSs 
of the (1, 0, 0) and (1,1,0) types remain stable, except in 
narrow regions close to edges of the corresponding band 
gaps. As a representative example, Fig. [3a) displays the 
evolution of the soliton of the (1, 0, 0) type, corresponding 
to point F9 in Fig. \Hb), which represents a GS built of 
1650 atoms with n — 10.3hoj±. While the latter is much 
greater than the transverse quantum hou±, the GS remains 
stable under the action of the 3D perturbations. The loosely 
bound GSs trapped in weaker OLs, such as the one in Fig. 
[2ti), are stable as well. 

Although vortices with the m — 2 are usually unstable 
against splitting 11711 , they may be stable in trapped config- 
urations Hi]. We have found that the GS of type (1,2,0) 
is stabilized by making it heavier: While it is unstable if 
built of < 300 atoms, the one shown in Fig. 0c) , which 



contains 450 atoms, is robust in the simulations. On the 
other hand, GSs with n r > 1 (that is, those featuring the 
complex radial structure) are unstable against quadrupole 
perturbations. In particular, Fig. Ob) displays the insta- 
bility for the GS of type (1, 0, 1), initiated by a very weak 
(0.1%) quadrupole deformation of the trapping potential, 
acting for 0.2 ms. GSs with n — 2 are also unstable, simi- 
lar to the previously mentioned "subfundamental solitons" 
in ID models [4]. Finally, the stability of GS complexes 
coincides with that of their fundamental constituents. 



IV. CONCLUSION 

Up to now, matter- wave GSs were actually created in 
shallow lattices (s ~ 0.7), with a weak radial confine- 
ment (E R /huj± ~ 44) [70. The respective linear spec- 
trum does not exhibit true band gaps. This can be checked 
using the approximate dispersion relation for the lowest- 
energy band in the shallow OL, E(q)/E R — (qd/n — 

l) 2 - V 4 MA - 1 ) 2 + 5 7 16 El- As follows from this, 
the width of the lowest band, E(l) - E(0) ~ 0.83E R , is 
much larger than the gap, huj± ~ 0.027^, which sep- 
arates the replicas corresponding to higher values of az- 
imuthal number m, and hence different Bloch bands over- 
lap in this case. (For parameters of the experiment from 
Ref. [7], true gaps open up only at s > 16.) Therefore, 
the solitary modes obtained under such conditions are ac- 
tually quasisolitons, decaying on a time scale of ~ 60 ms 
|7|]. The present results suggest possibilities for the cre- 
ation of robust fundamental and compound GSs, including 
vortical ones, which place themselves in true band gaps. 
Simultaneously, this approach opens possibilities for cre- 
ating truly multidimensional matter-wave solitons, which 
thus far have been elusive in experiments. 

In addition to the BEC, the proposed scheme may be 
used for the creation of 3D "light bullets" of the GS type in 
self-defocusing optical media with the anomalous group- 
velocity dispersion. In that case, t and z in Eq. (Q]) play the 
roles of the transmission distance and reduced time 111311 . 
The corresponding lattice cannot be created as a material 
structure, but it may be induced by a beating wave launched 
at a different wavelength, similar to the virtual grating used 
for the creation of optical GSs in ID [19]. Thus far, no 3D 
solitons have been created in optical media. 
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